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Abstract 

The high temperature phase of the electroweak standard theory is described by a strongly 
coupled S'C/(2)-Higgs-model in three dimensions. As in the Abbott-Farhi-model Higgs and 
M^-boson are low lying bound states. Using a method by Simonov based on the Feynman- 
Schwinger representation of correlators we calculate the masses of these states. Our results 
are compared with lattice masses. 
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The electroweak standard theory successfuhy predicted massive vector bosons making 
use of the Higgs mechanism. Much of its strength is based on the possibihty to do rigor- 
ous perturbative calculations. Later on interest also turned towards the behavior of this 
theory at high temperature that has been realized in the early universe. The proposal 
that there exists a hot phase where the Higgs field obtains a positive thermal mass and 
the SC/(2)i(gi?7(l)-symmetry is restored [Q] is commonly accepted. Investigations of the 
electroweak phase transition have to treat the hot and the cold phase at the transition 
temperature in an appropriate way. 

Perturbation theory in the broken phase has been studied extensively For 
A(/?^-couplings corresponding to zero temperature Higgs masses below 70 GeV it works 
very well. A perturbative treatment of the hot phase would predict a vanishing Higgs 
vacuum expectation value and a massless gauge boson. However perturbation theory 
alone is not a good guide in this phase due to IR-problems. 

The electroweak standard theory at the transition temperature allows the high tem- 
perature expansion. To good accuracy it can be described by an effective 3-dimensional 
5C/(2)-Higgs-model whose parameters can be calculated from the fundamental 4-dimen- 
sional theory All IR-problems appearing in the unbroken phase are represented by 
this effective theory. It might be compared to 3-dimensional QCD but with a fundamen- 
tal scalar field $ instead of quarks and a 5C/(2)-gauge group. $ corresponds to the 
perturbative Higgs particle. We avoid this name to preclude confusions. Its mass will 
be called Lagrangian mass. 

As in QCD the most secure way to calculate the properties of such a system is a Monte- 
Carlo simulation on the lattice. This is done intensely and first results are appearing 



|11|-|]1J]. Of course being aware of the long way to reliable hadron lattice results one has 
to question if the present lattice sizes are appropriate for measuring certain quantities. 
An interesting alternative way to such 3-dimensional systems at high temperature is the 
exact renormalization group approach It requires, however, simplifying assumptions 
for practical use, and these are not easily controlled. Recently it was speculated from 
the investigation of 1-loop gap-equations in the two settings with small and large Higgs 
vacuum expectation value, that the hot phase is just another Higgs phase with different 



parameters [16|. But this does not agree well with the masses mw and mn observed on 
the 3-dimensional lattice [0, 14 1. The subject of this letter is to propose a model 
which explains these correlation masses; we compare with the lattice results at the end. 



With rising 3-dimensional gauge coupling in the infrared [17| it is reasonable to propose 
confinement like in QCD. As in the Abbott-Farhi-model [18| physical states are bound 
states and singlets with respect to the gauge group. They may nevertheless be in a 
nontrivial isospin representation. Let us start with the 5C/(2)-Higgs-model Lagrangian 

£ = TiFkiF'^' + pfc(/>)tpV) - ^('AV) • (1) 

Dk = dk — igAk is the covariant derivative with the gauge field Ak = (r* are 

the Pauli matrices), Fki = [Dk^Di] is the field strength tensor and (p = {^^) is the 
scalar isospin doublet. C is invariant under the gauge transformation U{x) G SU{2) 

Dk ^ UDkU^ <P ^ . (2) 

The scalar doublet can unambiguously be written as (p = ^[^), where the 2x2-matrix 
$ is a linear combination of l2x2 and r*, and ^Tr = (f^cj) . Using this notation 
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uncovers the SU{2)(><)SU{2)j -invariance of the Lagrangian (eq. (||)), namely $ — > U^V , 
of which the left SU{2) is gauged while the isospin SU{2)j is global. 

The confining theory has (among others) bound states corresponding to the local 
interpolating fields 

TV and Tr$^L>fc$r* . (3) 

The first one is a scalar isospin singlet and is identified with the "Higgs particle" . The sec- 
ond one is an vector isospin triplet, the 'W-boson" . The scalar triplet operator Tr $^<I> r* 
vanishes. The operator Tr ^^Df^^ does not correspond to a vector singlet bound state 
but is identical to dkTi ^"^^ and generates the scalar singlet at non-zero momentum. 
There are no bound states with these quantum numbers. 

The nonlocal versions of the operators in eq. (^) are proportional to 

Tr ^{xyT{x,x)^{x) and Tv <^{x)^T{x,x)<^{x) . (4) 

The space index of the vector operator is fixed by the direction of the link operator T (also 
known as gauge field transporter). The latter is defined as {V denotes the path ordering) 

T{x,x) =Vexp (^ig J Akdzk^ . (5) 

These are the direct continuum counterparts of the lattice observables interpreted as Higgs 
and VF-boson. 

It has been shown in the framework of sum rules, that this model reproduces the 
phenomenology of the standard theory in the Higgs phase without any use of naive spon- 
taneous symmetry breaking but with an appropriate pattern of Tr vacuum structure 



|19|-[^]. (There is no fundamental difference between the 3 and 4-dimensional case.) Of 



course already at that time there was the other option (the genuine Abbott-Farhi-model) 



that gluon condensates like in QCD lead to a strongly bound phase |2^. Indeed this is 
the picture we propose for the hot phase of electroweak theory. 

The emerging physical picture would thus be very similar to QCD in 3 space time 
dimensions, but with spinless constituents, the fundamental scalar fields. It is tempting 
to try to make connection to the confining quark model of QCD which in a very trans- 
parent way can explain many features of the bound state spectra. Since the fundamental 
scalar particles are light compared to the string tension a we cannot use the simple non- 
relativistic quark model. 



In the case of QCD Simonov |23] has used the Feynman-Schwinger representation of 
hadronic correlators in order to treat light quarks kept together by a linear confining 
potential. The method is easily adopted to the case treated here, i.e. a theory in 3 
space-time dimensions with scalar "quarks". In this case we even avoid the problems 
of the large spin-interactions which lead presumably to zero modes and chiral symmetry 
breaking. Another difference is that the PV^-boson is an isospin-triplet, a state which has no 
analogy in QCD. Simonov's method is nevertheless applicable since the isospin structure 
is not changed by the propagation. The propagator of the matrix field ^^{x)aa to ^{y)bi3 
is given in the Feynman-Schwinger representation by (a, /? are gauge and a, b are isospin 
indices) 

G{x,y)aabi3 = (isexp(-m|s) J Vz exp (^-^ zl{T)dT^ T{x,y)a(3 5ab ■ (6) 
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The path integral / Vz runs over all curves z(r) with tangent i(r) connecting x and y. 

Using a similar expression for the anti-scalar propagator from x to y one can put 
things together to get the Green function for a scalar-anti-scalar state at position (x, x) 
propagating to (y, y) in the quenched approximation 

G{x,x,y,y) = J ds j ds j j PzPzexp -|- s)^ x (7) 

X exp (^-^ zl{T)dT - ^ fkiT)df^ V exp (^ig (j) _ Akdzk^ . 

The last factor of eq. ^ is just the Wegner Wilson integral over the closed loop formed by 
the gauge field transporters T of the propagators (eq. (^)) and the bound state operators 
(eq. (^)). In a first approach it is approximated by the area law 



V exp ( ig (p A^dz^ ) oc exp {—aF) , 

\ J x,x,y,y / 



(8) 



where F is the area and a the string tension. 

We are evaluating the Green function eq. (0) following Simonov |23|. Making some 



plausible assumptions on the (nearly) minimal surfaces the surface F is parameterized 
by Zk{T) and Zfc(r). In a next step the "center of mass trajectory" Rj. and the relative 
coordinate Uk are introduced 

ss fl l_\ 

Rk = — — -Zk + -Zk Uk = Zk- Zk . (9) 



s + s \s s 

The path integrals over z^ and z^ transform into path integrals over and Uk- 

Making the assumption that the center of mass motion is dominated by the classical 
path the path integral over R^ can be replaced by the parameterization Rkij) ■ Without 
loss of generality we choose Ri = R2 = , R3 = 7© = ^ with < 7 < 1 , where 
is the distance between the centers of mass of the bounded scalar-anti-scalar system 
at position (x,x) and at position {y,y)- This assumption is the most stringent one and 
should become better for heavier states (as compared to \/o'). 

One then can rewrite the Green function (eq. (0)) as (for detail see ref. ||23t| ) 

G{x,x,y,y) oc / ds ds Vuexp{—B) (10) 
JO JO J 

with the three dimensional action 



re 

B = I d^ 







^ e / 4ss A{s + s)\d^J V i ^ 



ill) 



It is convenient to introduce the new variables 

2s 2s Hi + fi2 

The ui and U2 path integral can be substituted by the solutions of the two dimensional 
Schrodinger equation 

L ^ 

2fl \dui 



Hil){ui,U2) = e'4){ui,U2) with H = - 7^ ( ^ + ^ ) + o-y-uf + ^2 . (13) 
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a-ni 


1=0 1=1 1=2 


n=l 
n=2 


1.74 2.87 3.82 
3.67 4.49 5.26 



Table 1: The dimensionless eigenvalues o„i of eq. (|15| 



This is most easily seen if one interprets B (eq. ([Tl|)) as the the Euclidian action of a point 
particle in two space and one time dimension. 

The last path integral over ii3 is now trivial and gives a negligible contribution. 

The remaining s and s integrals are finally evaluated by the method of steepest descend. 
The correlation mass M of the bound states, defined by G oc exp(— MG) , is hence found 
by minimizing the integrand of B evaluated at the solutions of the Schrodinger equation 
with respect to and /i2- Since the Lagrangian masses of the both fundamental scalars 
are the same (ms) one gets only one condition 



9M(/x) 







+ ^^ + e{^l) 



(14) 



/i 

with fii = fi2 = fJ- = 2fl . e{fi) is the eigenvalue introduced in eq. (13). One may easily 
express e(/u) by the dimensionless eigenvalues a„/ of the equation 



- Hp)" 



aniHp) 



through 



enlip) 



a 



2/3 



P 



1/3 



(15) 



(16) 



Here n is the principal quantum number and / the orbital one. The numerical values of 
ttni for the lowest quantum numbers have been obtained numerically. They are listed in 
table 1. In order to compare with the potential and with the lattice size the argument of 



the wave functions has been scaled back in units {g^T) ^ 
(radial symmetric) solutions of the Schrodinger equation 

the potential in Fig. 



versus r 



From the steepest descend equation (U) we obtain 

li = ^fa z{anuml/af/'^ 
where z(a, y) is a solution of the cubic equation 



V'n=i,/=o and ijjn=i,i=i , the 
13|), are plotted together with 



(17) 



yielding 

z{a,y) 



1 

-a z 
3 



y 







(18) 



2i/3 



-1/3 

21 y + A/729y2 - Aa^ ] + 



3 21/3 



27 y + V 7292/2 - Aa^ 



1/3 



(19) 



The mass M„; of the bound state with the quantum numbers n and / is read off eq. ([T 
using eq. ([l6|) and eq. ([l 



2ml 



3/2 



(20) 
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ml /a 

Figure 1: The masses of the bound states vs. the squared Lagrangian mass. The full lines 
are the masses of the s-states, the dashed lines are those of the p-states and the dot-dashed 
line corresponds to the Id-state 

Mni of the lowest bound states is plotted in Figure |l] versus the mass of the fundamental 
scalar (which corresponds to the current quark mass in QCD). 

Our results show the expected hierarchy of states, i.e. the lowest lying state is an 
isoscalar S-state (the "Higgs"), followed by an isovector P-state (the "W-boson"). They 
should now be compared with existing lattice data . These data are given for certain 
values of lattice parameters f3c and Ph, where Ph essentially fixes the temperature, and 

the lattice constant is a = p^^g'ix ■ Masses in units of g^T are obtained by multiplying 

the lattice ones with 

We consider lattice data of ref. corresponding to a zero temperature Higgs mass 
of 35 GeV, and concentrate on the data point (5g = 12, f3H = 0.3412 where the static 
force has been measured as well. The mass values are mu = (0.65 it 0.1) g^T and 
m\Y = (1.75 ± 0.08) g'^T , measured on a 30^ lattice. The static force is consistent with 
a confining potential with string tension a = (0.13 it 0.02) (g'^T)'^ up to a distance of 
4:{g'^T)~^ where it has been measured. At some (large) distances screening should set in, 
which, however, is presumable not relevant for the lowest bound states. 

For a comparison, we have to estimate the corresponding Lagrangian mass ms of the 
fundamental scalar. From it is obtained as 

(rmy ^ [PgY / 2(l-2/?fl-3/3g) 

WTj U M Ph 
with Pfi = 'gi'^ ! aiid the subscript indicates that the expression has to be sub- 
tracted at the critical Ph , which is equal to 0.34138 in the present case. This leads 
to 7713 = O.lJg^T at Ph = 0.3412 . In eq. (|2l| ) we assume for simplicity that 
vanishes at the critical Ph- Strictly speaking, depends on an arbitrary normalization 
scale /U3; in deriving eq. (^) we adopt the natural choice fi^ oc g^T . 
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4.0 6.0 

r/(g2r)-i 

Figure 2: The potential and the bound states. The thick full line shows the complete 
potential while the dotted line is its linear contribution. The dashed lines are the energies 
e^i of eq. (^). They have been used as zeros for the bound state wave functions, which 
are in arbitrary units. We used a = 0.13 ((/^T)^ and = 0.17 g'^T for this plot. 



A comparison with lattice data can now be performed based on equations ( [20| ) (see 
also Fig. |lj). Using the string tension a = (0.13 it 0.02) [g^T)"^ and the Lagrangian 
mass m3 = O.lJg^T from the lattice data one sees that our model yields qualitatively 
the right results of a substantial splitting between these states. The predicted mass of the 
VF-boson {m\Y = (1.47 it 0.11) g'^T) agrees quite well with the lattice mass, in view of the 
approximations made. Taking into account finite-a effects the lattice value may come down 
|14| improving the agreement. The mass of the composite Higgs (m// = (1.06 it 0.07) g'^T) 
calculated by us is a factor 1.6 larger than the lattice value. 

One may also take into account the attractive Coulomb force. The force measured on 
the lattice is excellently fitted by 



0.123 (g'^T)'^ Kiil.Olg'^T r) + 0.13 (^^T)^ 



(22) 



(cf. Fig. 5 of ref. jljl). The potential is obtained from F by integration fixing it to ar at 
large distances. The first term corresponds to the exchange force of a particle with mass 
l.Ol^f^T , the constant part is the string tension. The influence of the Coulomb force is 
small: it lowers the mass of the composite Higgs by only 2%, that of the W by 0.2%. This 
is due to the fact, that the bound state solutions are much larger than the range of the 
Coulomb force (cf. Fig. The r-dependent part of the force can therefore be neglected 
in good accuracy justifying the use of the area law in eq. (|8|). 

As far as our model calculations are concerned various approximations are involved. 
Perhaps the most important one is neglecting the quantum fluctuations of the center- 
of-mass coordinates Rk. This approximation should work better for heavier (recurrent) 
bound states. The deviation of the Higgs correlation mass may be explained by this. 

The lattice results are influenced by finite volume as well as finite lattice spacing effects. 
Although finite size investigations did not show significant effects on the masses one should 
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keep in mind that the predicted size of the bound states is comparable to the lattice size, 
which in this case is 10{g^T)~^ with periodic boundary conditions (see Fig. |2|). 

Recurrent states in the isoscalar S-channel would be difficult to identify because they 
are predicted in the mass range of the VF-balls of pure SU{2) gauge theory (l.S^^T) |24|. 
More interesting are the d-states (1=2). The corresponding 1/V^-ball states have masses 
of the order of 2.5g^T which is significantly heavier than the predicted 1=2 scalar-anti- 
scalar bound state mass mid of about 1.8 g^T . Note, nevertheless, that this is above 
the threshold of the two Higgs channel corresponding to the measured masses. Lattice 
measurements of this channel are under way. 

The predicted masses are based on the evaluation of the Green function eq. (0), which 
could only be performed with some approximations. Even without these approximation 
one finds that the masses of the bound state model depend only on the string tension a 
and the Lagrangian mass m^. Considering lattice data in the whole Ph (i-e. m^) range one 
finds a stronger m3-dependence of m\Y/{g'^T) and mH/{g'^T) than we would predict if 
^/ajicp'T) was constant. This might indicate that y/a/{g'^T) itself depends significantly 
on ms, dropping as approaches zero. It would be interesting to verify or disprove this 
by lattice studies. y/a/{g'^T) could also show a sizable A-dependence (decreasing with 
increasing A). 

One should also keep in mind that the static potential will not continue to rise linearly 
at larger distances but charge screening will set in. The screening length has not been 
measured so far. As a consequence of charge screening, heavier bound states will disappear. 
The mass of a bound state close to the threshold would be affected as well. 



In conclusion the hot electroweak phase appears as an intriguing case of a confining 
3-dimensional gauge system. It is similar to 3-dimensional QCD but with a fundamental 
scalar $ instead of quarks and hence without spin-interactions and without the problems of 
chiral symmetry breaking. The masses of the low lying bound states have been calculated 
from the Green function within some approximations. The comparison with lattice results 
is satisfactory. Our discussion is open for refinement both on the lattice side because 
of the limited lattice size and on the side of calculating the masses of relativistic bound 
states, a notorious problem also in QCD. For the latter question also other methods can 
be envisaged, including sum rules for the 3-dimensional Abbott-Farhi-model [0]-|21]. An 
enlarged version of this letter is in progress p5||. 
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